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ATOMIC  STRUCTURE  IN  STRONGLY  COUPLED  NEON  PLASMAS 


I.  Introduction 

Since  it  is  now  possible  to  produce  very  hot  plasmas  at  greater  than 
solid  density  in  the  laboratory,  it  is  meaningful  to  construct  theoretical 
models  of  such  systems  that  allow  for  accurate  determination  of  the 
systems'  properties.  For  calculations  of  atomic  properties  the  usual 
approach  is  to  iteratively  solve  a  set  of  coupled  equations  statistically 
describing  the  charge  distributions  and  an  effective  electron-ion 
interaction  potential.  Incorporating  this  potential,  the  bound  and  free 
electron  distributions  are  found  from  the  Schrodinger  equation.  Thus, 
given  an  ion  where  the  bound  orbits  are  not  externally  specified,  the 
solution  of  the  equations  directly  gives  orbital  energy  eigenvalues  and 
(fractional)  populations.  The  wave  functions  and  the  effective 
electrostatic  potential  obtained  in  this  manner  can  be  used  to  find 
spontaneous  decay  rates  and  cross  sections  for  various  atomic  processes 
characterizing  radiation,  the  spectrum  of  which  can  be  employed  to  diagnose 
the  plasma  environment. 

Thomas-Fermi  and  Hartree-Fock  statistical  models  have  been  applied  to 

highly  ionized  atoms  in  dense  plasmas  and  subsequently  applied  to  a 

strongly  coupled  neon  plasma. However,  ion  correlations  were  neglected 

•• 

in  these  approaches.  A  self-consistent  set  of  Schrodinger-  Poisson 
equations  including  ion  correlations  was  developed  by  Skupsky^^  to  study 
the  plasma  microfield  effects  on  a  high-Z  impurity  ion  embedded  by  a  dense 
fully  ionized  low-Z  plasma.  An  improvement  over  this  method  -  the  quantum 
mechanical  treatment  of  the  free  electrons  -  was  made  by  Davis  and 
Blaha.(^)  In  a  similar  manner  density  functional  theory  (DFT)  has  been 
employed  to  investigate  level  shifts  and  screening  effects  in  the  impurity 
problem.  ^»8) 

The  Inclusion  of  ion  correlations  in  these  latter  models  is 
accomplished  using  a  Boltzmann  distribution  under  the  assumption  of  nearly- 
dasslcal  ion  Interactions.  In  the  case  of  the  one-component  plasma 
(dynamic  ions  in  a  neutralizing  background  charge),  the  assumption  of  a 
Boltzmann-like  form  for  the  ions  would  be  erroneous  for  values  of  the  ion 
coupling  parameter, 

_  (z)Vg 
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greater  than  about  three. Here  Z  is  the  effective  ionic  change,  rQ 
is  the  ion  sphere  radius,  and  3  *  l/k_T.  This  discrepancy  is  not  as 

O 

significant  for  a  "real”  two-species  plasma  because  the  mobile  electron 
fluid  is  able  to  provide  more  effective  screening,  but  has  yet  to  be 
investigated  in  the  two-species  model  for  T  >  2  and  for  ions  other  than 
hydrogen.  If  one  can  utilize  a  model  that  is  expected  to  provide  accurate 
distributions  for  a  strongly  coupled  system,  one  can  also  use  that  model  to 
examine  the  validity  of  using  self-consistent  statistical  models  in  the 
strongly  coupled  regime. 

Implicit  in  all  the  methods  discussed  here  is  the  assumption  that  the 

lifetime  of  the  ionic  state  is  long  enough  so  that  the  plasma  has  time  to 

be  polarized  by  the  ion.  Since  the  polarization  (correrlation)  time  is  of 

the  order  of  the  inverse  of  the  plasma  frequency,  all  of  the  cases  we  are 

considering  can  be  considered  long-lived  (a  typical  state  lifetime  —  the 

most  rapid  destruction  mechanism  being  collisional  de-excitation —  may  be 

of  the  order  of  10"^  -  lO-*^  s;  u  *  is  about  10-^  g).  Each  model  also 

P 

assumes  that  since  the  ion  state  exists  through  many  plasma  periods,  the 
concept  of  a  time-averaged  potential  for  atomic  calculations  is  meaningful. 

It  should  be  noted  that  the  calculations  involving  the  self-consistent 
method  presented  below,  like  those  of  Ref.  8,  form  an  "average  atom"  ion 
model,  in  contrast  to  models  which  detail  the  ionic  configuration.  The 
eigenenergles  obtained  below  do  not  represent  the  spectrum  of  an  ion  in  a 
specific  configuration  (i.e.,  hydrogenic),  but  represent  an  average  over 
many  ions  in  partially  ionized  states. 

We  investigate  here  the  energy  eigenvalues,  charge  distributions,  and 

effective  electron-ion  potentials  for  strongly  coupled  neon  plasmas  using  a 

self-consistent  DFT  model  similar  to  that  described  in  Ref.  8.  These 

results  are  compared  with  those  obtained  from  the  solution  of  the  two- 

component  plasma  hypernetted  chain  (HNC)  equations,  which  are  assumed  to  be 

valid  at  these  densities  and  temperatures.  The  results  will  indicate  the 

•• 

inadequacy  of  the  Debye-Huckel  (DH)  and  ion-sphere  (IS)  calculations  when 
13  i.  T  <  1  . 

In  Sec.  II  we  describe  the  atomic  and  plasma  models  we  will 
consider.  Section  III  contains  the  results  of  computations  employing  these 
models.  The  results  are  discussed  in  the  concluding  section. 

2 


i 


* 


II .  Models 


We  consider  an  ion  of  nuclear  charge  Z  in  a  plasma  in  which  the 
average  effective  charge  is  Z.  Z  is  equal  to  Z  minus  the  mean  number  of 
bound  electrons  per  ion  and  is  a  result  of  the  model.  Density  functional 
theory  leads  to  a  system  of  equations  that  must  be  solved  self- 
consistently.  The  electrostatic  potential  is  given  by  the  Poisson 
equation. 


V(r)  - 


Ze  ,  ,  2  f  1  /  r,  2,  .  ,  N 

T"  +  4ire  1  7  to  dr  r  <pe4Vpi) 

+  /*  dr  r  (pe+pb+p1)  ]. 


(1) 


The  plasma  is  assumed  to  be  in  thermal  equilibrium  and  all  electrical 
charge  distributions  are  assumed  to  be  spherically  symmetric.  In  Eq. 
(1),  pb  is  the  local  density  of  bound  electrons 

pb  -  -  (4-rV1  I  *nl  P2l(r)  .  (2) 

n£ 


The  b^  are  the  state  occupation  numbers  (£  bn^  ■  Z  -  Z  )  and  Pn^(r) 

are  the  radial  wave  functions  found  from  solving  the  Schrodinger  equation 

where  the  Interaction  potential  is  V(r)  from  Eq.  (1)  with  p.  set  equal  to 

D 

zero, 


£(*+1) 


+  V(r)  -  Vxc(r)  +  Vxc(»)  +  En,]  P^r)  -  0  , 


(3) 


where  Ent  is  the  energy  eigenvalue  of  state  n£.  Here, the  exchange- 
correlation  potential,  ^xc(r),  has  been  calculated  by  Gupta  and 
Rajagopal^1^  and  by  Dharma-wardana  and  Taylor. 

p*  is  the  local  charge  density  of  free  electrons.  It  is  represented 
by  a  Fermi-Dirac  energy  distribution  beyond  a  spherical  boundary  large 
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enough  so  that  the  plasma  at  the  boundary  may  be  considered  neutral. 

Inside  this  sphere  the  free  electrons  may  be  treated  quantum-mechanically 

and  are  described  by  wave  functions  that  are  solutions  of  the  time- 
«• 

independent  Schrodinger  equation,  i.e. 


where 


P  -  -  /"dr  W(k)  -i-y  l  (24+1)  Fj*  (r) 
6  °  k  r  £-0  ™ 

2  2 

W(k)  »  ^  [1  +  exp  {(^2  “  U)  3}]"1; 
ir 


(4) 


(5) 


here  y  is  the  chemical  potential  of  the  free  electron  gas  determined  from 


J"  dk  W(k)  -  p  («)  -  -  n  , 
o  e  e 


(6) 


where  n@  Is  the  mean  electron  density.  The  free  electron  wave  functions, 
F^  are  solutions  of  Eq.  (3)  with  the  replacement  of  the  eigenenergy,  En^, 

O 

by  the  electron  kinetic  energy,  k  .  The  ion  charge  density  is  assumed  to 
take  the  Boltzmann  form 


pi 


\  e~SV(r) 

Z 


(7) 


At  r  ■  pj  ■  -  p^,  insuring  neutrality;  we  also  have  the  boundary 

conditions  rV(r)  ♦  0  and  0  as  r  approaches  infinity.  Equations  (1-4) 

and  (6-7)  are  solved  self-consistently  with  these  boundary  conditions  to 
yield  En£,  Pe>  pt,  and  V(r). 

In  order  to  gauge  the  reliability  of  the  above  model  in  a  strongly 

coupled  plasma  we  turn  to  a  semiclasslcal  treatment  of  particle 

correlations  that  has  been  found  to  accurately  reproduce  molecular-dynamics 

calculations  in  this  regime.  In  this  approach  -  the  two-component  plasma 

(TCP)  -  the  ions  and  electrons  are  treated  as  classical  particles  that 

Interact  through  effective  two-body  potentials  which  deviate  from  pure 

Coulomb  behavior  at  short  distances  such  that  the  essential  quantum 

diffraction  effects  are  simulated.  A  particular  form  has  been  suggested  by 

Deutsch^*^  and  used  in  the  computer  simulations^^ .  This  form  uses  the 

reduced  mass  de  Broglie  wavelength,  X  . ,  where  a  and  0  are  species  labels, 

up 

as  a  quantum  mechanical  cutoff  parameter,  i.e. 


(8) 


Wr) 


?a?6e 


[1  -  exp  (-r/^ag)]; 


is  the  charge  of  species  a  and 


«/(2,Ua6k8T>‘ 


where  u  „  is  the 
ocB 


reduced  mass.  This  potential  is  finite  at  the  origin  and  is  expected  to 
give  reasonable  results  for  nondegenerate  plasmas  so  long  as 


X  /r  «  1  (  X  is  the  smallest  of  the  three  X  „  .)  This  condition  is 

ee  o  ee  _  „  aB 

equivalent  to  T  «  9  (Z;  /(T^)1  . 

In  order  to  include  the  plasma  many-body  effects,  the  binary  inter¬ 
actions  defined  in  Eq.  (8)  are  used  in  the  hypernetted  chain  (HNC) 

equations.  This  is  an  approximate  integral  equation  method  for 

calculating  static  correlation  functions  for  systems  of  particles  with  long 
range  potentials  and  has  proven  to  be  accurate  for  strongly  coupled 
hydrogen  plasmas. The  quantities  of  interest  are  the  radial 


distribution  functions  (rdf's),  g  „(r), 

Up 


which  contain  the  static 


structural  information  in  the  TCP.  The  HNC  approximation  for  the  rdf's  is 


gaB(T)  -  exp  [-8VaS(r)  +  ha6(r)  -  c^r)]. 


where  the  total  correlations 


ha6(r)  "  8a8(r)  "  1 


are  related  to  the  direct  correlations  c^g  by  the  Ornstein-Zernicke 


equations 


ha8(k)  "  Ca0(k)  +  ^  hcrr(k)  cYB(k)* 


Here  the  Fourier  transform  is  defined  as 


h  ft(k)  "  fn  dr  r 
aB  a  '  o 


2  slnkr 


Wr)- 


Eqs.  (9-12)  are  solved  iteratively  for  o,B  ■  i,e.  The  rdf's 

•• 

generated  by  this  procedure  reduce  to  their  Debye-Huckel  (DH)  forms  in  the 


limit  of  weak  coupling  (T  «  1),  but  are  considerably  different  from  the 
DH  approximation  when  r  is  order  one  or  larger. 

The  TCP  is  a  model  system  of  point  charges,  ions  with  charge  +  Z  and 
free  electrons  with  charge  -1.  Formally  the  HNC  scheme  requires  the  exact 
Zf  as  an  input  parameter;  this  is  necessary  if  the  ionic  and  electronic 
distribuiton  functions  are  to  be  examined.  In  order  to  find  the  effective 
potential,  however,  only  a  rough  guess  of  Z  will  suffice  to  determine  much 
of  the  V(r)  curve. 

The  effective  electron-ion  potential  and  the  screening  function, 
e  (k),  are  defined  via  Poisson's  equation,  the  Fourier  transform  of  which 
is  given  by 


V(k) 

HNC 


4tre2Z 

k2e(k) 


[~siLW  -  sle(k)//  z  ] 


(13) 


The  static  structure  factors  are  defined  by 


Vk>  -  5.6  +  (lta:8l>1/2  Vk>  ‘ 


(14) 


Close  to  a  test  point  ion  of  charge  Z  the  free  electron  distribution 

determines  V(r)  ;  the  ion-ion  rdf  is  negligible  out  to  a  distance  of 
HNC 

about  one-half  rQ.  In  this  region  Poisson's  equation  is 


vV(r)  -  4ire2  Z  [«(r)  -  (r)],  r  £  r  /2,  (15) 

HNC  Z 

where  the  prime  on  V(r)  indicates  the  test  ion  has  charge  Z ,  not  Z.  For  a 

HNC  _ 

given  temperature  and  electron  density,  a  higher  value  of  Z  simply  pulls 

the  electron  distribution,  h^e(r),  in  tighter,  an  effect  that  essentially 

compensates  the  Z  prefactor  to  hi@.  The  result  is  that  the  function  in 

brackets  in  Eq.  (16)  is  very  nearly  Insensitive  to  the  mean  ionic  charge, 

1«  e. 


V'(r)  -  4ire2  Z  f(r)  (16) 

HNC 

where'  f(r)  is  a  function  nearly  indf  **endent  '  Z.  This  is  the  rationale 
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behind  the  form  of  the  potential  In  Eq.  (13);  V  (r)  is  a  screening 

HNC 

function  dependent  on  density  and  temperature  scaled  by  the  nuclear  charge 
Z.  We  find  that  this  form  very  nearly  reproduces  the  potentials  found  in 
the  quantum-mechanical  self-consistent  model  described  earlier. 

If  the  particle  distributions  are  required  (an  accurate  value  of  Z  is 
needed  to  obtain  the  actual  distributions)  two  steps  are  necessary.  First, 
a  guess  of  Z  is  made  and  the  HNC  equations  solved  for  S^^(k)  and  Sj_e(k)* 

A  "guess"  of  the  potential  is  then  found  from  Eq.  (13).  This  potential  can 
then  be  used  in  Eq.  (3)  to  find  wave  functions  for  all  bound  states.  The 
integrated  wave  functions  provide  a  new  Z,  in  which,  when  used  in  the  HNC 
code  a  second  time,  provide  a  new  potential  and  the  needed  distributions. 
Generally  only  the  one  such  iteration  is  required. 

The  definition^^  of  y(r)  in  Eq.  (13)  implies  a  form  of  the 
dielectric  function  significantly  different  from  that  obtained  using  the 
fluctuation-dissipation  theorm^2^  (FDT),  although  both  forms  reduce  to  DH 
forms  in  the  proper  limits.  For  dense  plasmas  VHNC  from  Eq.  (13)  agrees 
much  more  closely  with  results  from  the  ion  sphere  model  (described  below) 
and  with  Thomas-Fermi*  calculations,  as  well  as  the  "potential  of  mean 

Ze2 

force"  approximation,  V^r)  *  — —  In  [gei(r)J,  than  an  effective 

Interaction  derived  from  the  FDT.  In  fact  VFDT(r)  shows  screening  that 
everywhe*—  has  a  larger  magnitude  than  vDn(r)»  T^e  De^ye  potential  itself 
is  known  already  to  predict  excessive  screening  in  plasmas  where  the 
validity  of  the  DH  approximation  is  questionable.  A  plasma  in  a  near 
metallic  state  (where  the  ion  sphere  model  might  be  used)  shows  a  form  very 
similar  to  the  potential  defined  via  Poisson's  equation,  which  is 
qualitatively  and  quantitatively  distinct  from  an  effective  interaction 
derived  from  the  FDT. 

As  the  plasma  approaches  the  limit  of  a  solid  structure,  the 
ion-sphere  (IS)  approximation  becomes  more  valid.  The  ion-sphere 
model^*^)  assumes  complete  ion  shielding  within  an  ion-sphere  radius  by  a 
uniform  cloud  of  electrons.  Poisson's  equation  in  this  case  yields 

vis(r)  ■  2'2  I  7  -  -K-  <3  -  T-i>  1  •  (17) 

0 

For  extreme  densities,  VIS  should  be  approximately  valid. 
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III.  Results 


We  consider  a  strongly  coupled  neon  gas  plasma.  Table  I  summarizes 
the  conditions  under  which  the  runs  were  made,  the  value  of  Z  being  a 
result  of  the  self-consistent  (SC)  model.  All  cases  have  r  's  in  excess 
of  two.  We  note  that  both  the  HNC  and  SC  models  reduce  to  correct  Debye- 
Huckel  results  in  the  limit  of  weak  coupling  (T  «  1  ). 

The  ion  charge  density  from  the  self-consistent  model  normalized  to 
the  backround  density,  p^/pO*),  is  equivalent  to  the  ion-ion  radial  distri¬ 
bution  function,  gj^.  Fig*  1  displays  the  ion  distributions  resulting  from 
SC  solutions  for  the  T  m  2.2  and  T  *  4.9  cases.  These  figures  are  compared 
with  gj^  from  the  HNC  approximation  using  the  effective  binary  interaction 
in  Eq.  (8)  and  with  the  Debye  form 


DII,  . 

*n(t) 


-  2  _r^D 
exp  [  -Ze  e  ] 


(18) 


where 

*  4irnee^(Z+l)g  . 


(19) 


g^  shows  the  tendency  of  the  DH  approximation  to  excessively  screen 
the  ions  in  dense  plasma,  an  effect  previously  seen  in  the  OCP^9,1^  and 
the  TCP^^\  The  HNC  rdf  is  assumed  to  be  the  most  accurate  of  the  three 
representations,  because,  since  -  lO-'*,  the  ions  are  essentially 
classical  particles  and  the  computer  simulations  have  supported  the  use  of 
the  HNC  approximation  for  classical  systems.  In  spite  of  the  fact  that  p^ 
in  the  SC  method  -  Eq.  (8)  -  cannot  reproduce  the  oscillations  around  gj^  * 
1.0  for  2  <  r/aQ  <  4  in  the  larger  r  case,  the  agreement  between  SC  and  HNC 
even  at  V  “  4.9  is  very  good.  The  small  difference  between  these  two  forms 
is  not  expected  to  alter  the  effective  potential^®) ;  we  will  test  the 
significance  of  the  difference  below. 

Fig.  2  compares  the  electron  density  profile  (including  both  bound  and 
free  electrons)  provided  by  the  self-consistent  method  around  an  ion  with 
the  ion-electron  radial  distribution  function  produced  by  the  HNC  code  for 
T  ■  2.2.  The  profiles  are  very  close  for  r/aQ  >  0.25.  The  innermost  r- 
point  calculated  on  the  Fourier  transform  mesh  in  the  HNC  code  is 


r/aQ  *  0.125.  This  also  corresponds  to  the  innermost  r-mesh  point  of  the 
potential,  since  V(r*0)/e  ■  Z,  interpolation  between  r  *  0  and  the  first 

mesh  point  is  possible  for  ^hnC^1^*  Extrapolation  of  the  HNC  g^e  to 

smaller  radii,  however,  would  not  be  meaningful,  since  ^g/a  =  0.12; 

thus  quantum  mechanical  details  are  important  in  this  region.  Forrest 

Rogers  has  investigated  this  subject  for  hydrogen  and  few-times  ionized 

argon. Since  our  goal  is  a  many-body  effective  potential  with  which  to 
examine  average  atom  calculations,  we  find  that  the  present  model  is 
adequate. 

The  SC  effective  potential  is  a  consequence  of  the  solution  of  the 
model.  This  function  in  the  form  rV(r)  appears  in  Fig.  3  for  r  *  2.2  and 
F  *  4.9.  The  HNC/Poisson  potential  -  Eq.  (14)  -  is  also  presented.  The 
two  forms  are  seen  to  be  very  similar  in  both  cases  indicating  the  apparent 
validity  of  the  quantum  mechanical  model  even  at  very  high  densities.  The 
Debye  potential  reveals  much  stronger  screening  except  for  large  distances 
where  rV(r)  tends  to  zero  for  all  models.  The  ion-sphere 

approximation  is  included  for  comparison:  It  agrees  rather  well  with  SC 
and  HNC/Poisson  at  short  distances,  but  predicts  even  larger  than  Debye 
screening  farther  as  r  becomes  larger  -  a  tendency  very  much  distinct  from 
SC  and  HNC.  The  overall  form  of  the  IS  function  is  very  different  from  the 
exponential  behavior  of  the  DH,  SC,  and  HNC/Poisson  functions,  a  result  of 
its  constraint  of  fixed  ionic  volume. 

Having  now  seen  that  the  self-consistent  formalism  can  provide 

reasonable  results  (compared  with  the  HNC  data)  for  these  strongly  coupled 

plasmas,  we  now  look  at  the  energy  levels  of  the  neon  ions.  Table  II  is  a 

compilation  of  negative  energy  eigenvalues  arising  from  the  solution  of  the 
•* 

Schrodlnger  equation  -  Eq.  (3)  -  within  the  method.  All  negative  (bound) 
energies  are  noted.  The  less  deeply  bound  or  absent  DH  values  (resulting 
from  more  severe  screening)  as  well  as  eigenvalues  found  by  using  vIS(r) 
and  V|}jig(r)  are  presented  for  comparison. 

As  a  test  of  the  significance  of  the  difference  between  the  two  forms 
of  the  ion  distribution  functions  -  the  SC  [  Eq.  (7)  ]  and  HNC  [  Eq.  (9)  ] 
a  run  of  the  SC  model  was  repeated  for  r  ■  3.4  using  g^C  as  a  fixed 
function  instead  of  Eq.  (7).  Those  figures  are  set  in  parentheses  in  Table 
II.  The  difference  is  indeed  minor  and  of  the  order  of  the  numerical 
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accuracy  of  the  coded  formalism. 

As  an  example  of  a  neon  plasma  at  extreme  conditions,  we  examined  the 
case  in  which  ng  ■  5»  lO^cm  ^  and  T  *  210  eV,  giving  a  T  of  13.1.  In  this 
regime  one  expects  to  see  considerable  difference  between  the  profiles 
produced  by  the  SC  and  HNC  methods.  In  Fig.  4  the  ion  distributions  of  the 
HNC,  SC,  and  DH  theories  are  reproduced.  The  HNC  rdf  shows  that  in  this 
case  ion  correlations  are  not  Boltzmann-like.  The  non-negligible 
oscillation  about  g^  -  1  shows  that  there  is  now  a  strong  indication  of 
ion  ordering.  Since  the  SC  method  utilizes  Eq.  (7),  this  effect  is  not 
seen  in  those  results.  Quantitatively  the  SC  and  HNC  rdf's  are  more 
dissimilar  than  the  less  coupled  cases,  although  the  DH  curve  is 
considerably  more  distinct  from  both  of  these. 

Fig.  4  indicates  that  the  effects  of  non-Boltzmann  ion  correlations 
are  expected  to  be  seen  in  the  potential  only  at  distances  of  r/aQ  greater 
than  one.  Inside  this  radius  the  SC  and  HNC  ion  distributions  are  similar 
enough  that  the  potential,  which  here  depends  on  the  electron  distribution, 
is  not  expected  to  be  greatly  affected.  The  Debye  potential  is  expected  to 
be  overly  screened  again.  Fig.  5  provides  the  calculated  potentials 
for  T  »  13.1.  The  ordering  seen  ir.  Fig.  4  is  obviously  manifested  as  the 
r-space  oscillations  in  V(r).  For  purposes  of  atomic  calculations,  this 
effect  will  have  little  significance  as  the  spatial  extent  of  the  Is  wave 
function  is  limited  to  the  volume  Inside  r/aQ  -  0.15.  For  comparison  we 
have  plotted  the  ion  sphere  potential  (crosses),  which  is  seen  to  coincide 
closely  with  the  SC  and  HNC  effective  potentials  up  to  r/aQ  “  0.5. 


IV.  Discussion 


Our  primary  goal  is  to  investigate  the  applicability  of  a  self- 
consistent  (SC)  quantum  statistical  method  to  strongly  coupled  neon 
plasmas.  As  points  of  comparison  we  Include  potential  calculations  from 
Debye-Huckel  (DH-— correct  for  T  «  1)  and  ion-sphere  (IS — assumed  correct 
for  T  »  1)  approximations.  The  solution  to  the  hypernetted  chain  (HNC) 
equations  incorporating  a  semiclasslcal  binary  pseudopotential, which 
has  been  found  to  be  accurate  in  strongly  coupled  hydrogen  plasmas, is 
the  plasma  model  whose  statistical  properties  the  SC  method  must  mirror  in 
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order  to  be  considered  valid  in  this  regime. 

The  SC  method  incorporates  ion  correlations  via  a  Boltzmann  factor 
with  the  self-consistent  potential  in  the  exponent.  Although  this  form  is 
approximate  and  cannot  predict  strong  correlations  which  result  in  spatial 
oscillations  in  gj^»  the  SC  ion  distribution  is  very  close  to  the  HNC 
profile  for  all  cases  considered  with  the  exception  of  T  =  13.  The  DH 
profiles  predict  more  closely  packed  ions  due  to  considerably  more 

screening,  a  charateristic  of  DH  theory  outside  of  its  range  of  validity 
and  known  to  be  incorrect. 

The  IS  profile  (not  depicted  on  Figs.  1  and  4)  is  a  step  function  at  r 
■  rQ  with  amplitude  zero  inside  this  radius  and  amplitude  one  beyond  rQ. 
This  extreme  form  is  not  appropriate  for  the  lower  T  cases,  but  is  nearly 
correct  for  T  ■  13  (rQ«0.64ao).  Of  course  the  structural  oscillations  are 
absent  in  the  IS  model  but  this  should  not  be  very  important  for 
calculations  involving  bound  electrons. 

The  inner  region  of  the  effective  potential  is  determined  mainly  by 
the  electron  distribution  around  the  ion.  The  HNC  and  SC  electron 
distributions  are  seen  in  Fig.  2  to  be  close  except  very  near  the  ion  where 
the  HNC  solutions  cannot  be  found.  The  ion  sphere  profile  here  possesses 
no  structure,  simply  a  horizontal  line  at  gle  **  1;  this  difference  is 
crucial  when  developing  the  effective  potential  that  is  used  to  investigate 
the  atomic  structure  of  the  ion. 

The  effective  potential  is  a  consequence  of  the  solution  of  the  SC 
method.  The  HNC  interaction  is  derived  from  a  screening  function  that  is 
nearly  Independent  of  the  average  effective  charge,  Z  (for  a  given 
temperature  and  density),  scaled  to  the  nuclear  charge  of  the  ions.  The 
derivation  of  vmjC^r^  ^rom  Poisson's  equation,  not  the  fluctuation 
dissipation  theorm  as  explained  in  the  last  section. 

The  SC  potential  is  found  to  be  very  similar  to  ^^(r)  for  the  lower 
T  cases.  The  HNC  potential  is  less  screened  than  the  SC  potential  in  the 
region  r  i.  rQ/2,  where  the  electron  distribution  essentially  determines  the 
form  of  the  potential.  The  electron  ”pile-up”  near  the  nucleus  is  larger 
in  the  SC  model,  more  effectively  screening  the  positive  charge.  The  SC 
electron  distribution  and  thus  the  effective  potential  in  this  region  are 
probably  more  accurate  than  the  HNC  results.  For  larger  r,  however,  the 
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ion  distribution  begins  to  effect  the  potential.  The  ion-ion  rdf  curves  in 
Figs.  1  and  4  show  the  ions  generally  less  packed  in  the  HNC  approximation 

than  in  the  SC  method,  evidence  of  the  greater  "pile-up",  but  the  structure 

is  not  simple.  The  enhanced  (non-Boltzmann-like)  ion  correlations  shorten 
the  range  of  the  calculated  potential.  The  HNC  potential  in  this  region  is 
probably  the  more  accurate  of  the  two. 

Fig.  5  indicates  the  presence  of  very  strong  enhanced  correlations 
effecting  the  potential.  Using  the  HNC  ion  distributions  in  the  SC  method 
in  place  of  the  Boltzmann  from  does  not  allow  g|^(r)  to  readjust  to  changes 
in  the  electron  distribution  in  each  iteration  (unless  the  HNC  code  is 

coupled  directly  to  the  SC  scheme,  a  procedure  we  have  not  undertaken). 

Only  in  the  T  •  13  case  might  this  be  important.  A  potential  constructed 
of  the  effective  potential  from  the  SC  theory  inside  r  %  r  /2  and  Vuvr-(r) 
outside  this  point  should  suffice  for  most  atomic  calculations.  For  bound 
state  energies  at  r  ■  13,  the  differences  between  HNC  an  SC  are  not 

important  (see  Table  II).  Calculations  of  continuum  properties  may  be 

effected.  This  problem  will  be  investigated  in  the  future. 

The  inner  region  is  more  important  for  obtaining  information  on  bound 
states.  Here,  the  ion  sphere  approximation  is  hampered  by  the  assumption 
of  uniform  electron  density  with  the  result  the  VIS(r)  predicts  more  deeply 
bound  states  than  HNC  or  SC  (which  are  in  turn  more  deeply  bound  than  the 
Inappropriate  DH  values).  This  is  true  even  at  the  extreme  case  of  T  * 
13.1  as  indicated  by  Table  II.  At  the  temperatures  and  densities  for  the 
cases  listed,  the  HNC  and  SC  potentials  predict  the  same  number  of  bound 
states  with  approximately  the  same  energies.  The  IS  approximation  predicts 
more  deeply  bound  inner  levels  but  may  result  in  less  energetically  bound 
outer  levels  (as  in  T  ■  2.2)  if  the  wave  functions  extend  into  the  region 
where  the  range  of  the  IS  potential  is  foreshortened  by  its  definition  of  a 
fixed  ionic  volume.  The  Debye  values  all  indicate  more  shallow  states, 
because  of  large  screening.  (If  in  the  definition  of  X^  in  Eq.  (19),  Z 
were  used  instead  of  Z,  the  result  would  be  even  more  severe  screening. 
Setting  Z  ■  0  in  Eq.  (19)  i.e.,  using  the  electron  Debye  length,  produces  a 
potential  devoid  of  any  ionic  contributions  to  the  correlation  functions. 
The  result  here  is  a  potential  that  lies  much  above  all  of  the  curves;  this 
approximation  provides  too  little  screening.) 


Since  the  HNC/Poisson  potential  appears  to  be  accurate  for  atomic 
calculations  in  neon  for  F  ^  2.2  (by  Fig.  3  and  Table  II),  this 
approximation,  which  is  very  easy  to  generate,  can  be  used  to  examine  other 
properties  of  such  strongly  coupled  systems  requiring  a  many-body 
potential.  For  calculations  at  higher  coupling  and  details  of  electron 
distributions  very  close  to  the  nucleus,  the  self-consistent  method  is 
needed. 

We  have  shown  the  self-consistent  model  produces  reliable  results  for 
strongly  coupled  plasmas  compared  to  hypernetted  chain  results  in  neon  up 
to  T  of  order  five  and  higher  if  HNC  ion  distributions  are  employed.  In 
addition,  we  have  shown  that  the  HNC  method  of  generating  correlation 
functions  provides  an  effective  potential  that  can  be  used  in  calculations 
of  atomic  properties  up  to  V  of  order  two  (for  neon).  Debye-Huckel  theory 
is  not  a  meaningful  approximation  in  strongly  coupled  plasmas.  Nor  can  we 
recommend  the  use  of  the  ion  sphere  potential  for  any  of  the  cases  examined 
here. 
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Table  I 


I.  Summary  of  selected  neon  plasma  conditions  described  by  the  models, 
is  the  mean  charge  per  ion.  r  is  the  ion  coupling  parameter. 


Electron  Density 

Temperature 

Z 

(cm-3) 

(eV) 

CM 

O 

H 

400 

8.77 

2.1024 

250 

7.61 

5.1024 

250 

7.85 

5.1025 

210 

8.02 

t  ,  \  ■>  t.  . 


•  •  * 


<»  « 
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TABLE  II 

Energy  eigenvalues  in  atomic  units  of  neon  plasmas  at  T  -  2.2,  3.4, 

HNC  SC 

4.9,  and  13.1  from  self-consistent  model  (using  in  place  of 

for  the  parenthetical  values  under  T  -3.4),  and  the  HNC/Poisson  model, 
Debye  model,  and  ion-sphere  model  potentials.  All  bound  level 
energies  are  given. 


Is  -42.3  -43.1  -43.8  -39.7 

2s  -6.15  -6.35  -6.45  -4.40 

2p  -5.93  -6.24  -6.40  -4.00 

3s  -0.64  -0.57  -0.23  -0.16 

3p  -0.50  -0.44  -0.12  -0.01 

3d  -0.22  -0.18 


-34.6  (-34.6)  -39.4  -41.8  -33.7 

-2.88(-2.89)  -3.77  -4.65  -1.54 

-2.31(-2.41)  -3.41  -4.55  -0.77 


-  13.1 


IS 

DH 

-27.1 

-0.66 

Figure  1 

Ion-Ion  radial  distribution  functions  for  (a)  r  ■  2.2  and  (b)  r  -  4.9 
in  HNC  (dots),  SC  (solid  line),  and  DH  (dashed  line)  approximations. 
The  distance  r  is  in  units  of  the  Bohr  radius,  aQ. 


r/a0 

Figure  2 


Electron  distribution  around  an  ion  for  r  ■  2.2.  Solid  line  is  the 
SC  model;  dashed  is  the  HNC  approximation.  The  smallest  radius  used 
in  the  Fourier  transform  within  the  HNC  scheme  was  r/aQ  ■  0.125. 
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Figure  3 

Effective  electron-ion  potentials  resulting  from  the  SC  model  (solid), 
HNC/Poisson  model  (dots),  ion-sphere  model  (crosses),  and  DH  theory 
(dashes)  for  (a)  F  *2.2  and  (b)  r  •  4.9. 
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Figure  5 

Effective  electron-ion  potential  for  V 
identical  to  Fig.  3. 


13.1  neon.  The  legend  is 
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